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ABSTRACT 

The radio-quiet quasar SDSS J1536+0441A shows two broad-line emission systems, recently inter- 
preted as a binary black hole (BBH) system with a subparsec separation; as a double-peaked emitter; 
or as both types of systems. The NRAO VLBA was used to search for 8.4 GHz emission from SDSS 
J1536+0441A, focusing on the optical localization region for the broad-line emission, of area 5400 mas 2 
(0.15 kpc 2 ). One source was detected, with a diameter of less than 1.63 mas (8.5 pc) and a brightness 
temperature Tj, > 1.2 x 10 7 K. New NRAO VLA photometry at 22.5 GHz, and earlier photometry at 
8.5 GHz, gives a rising spectral slope of a = 0.35 ± 0.08. The slope implies an optically thick syn- 
chrotron source, with a radius of about 0.04 pc, and thus ~ 5 x 10 10 K. The implied radio-sphere 
at rest frame 31.2 GHz has a radius of 800 gravitational radii, just below the size of the broad line 
region in this object. Observations at higher frequencies can probe whether or not the radio-sphere is 
as compact as expected from the coronal framework for the radio emission of radio-quiet quasars. 

Subject headings: black hole physics — quasars: individual (SDSS J153636. 22+044127.0) — radio 
continuum: general 



1. MOTIVATION 

Binary black hole (BBH) systems with subparsec 
scales are predicted in merging scenarios for galaxy evo- 
lution and also factor prominently in predictions for 
the g ravitational wave background (e.g., IColpi fc Dottil 
2009). But can such BBH systems be found? 
A promising search method - identifying candidate 
BBH systems through their optical broad-line emis- 
sion properties - recently yielded three such candi- 
dates: SDSS 
2009), SDSS 



J153636.22+044127.0 (Boroso n fc Lauer 
J0 92712 .65+294344.0 (jBogdanovic et al 



20091 iDotti eFani2"009h . and SDSS J105041.35+345631.3 
(Shi elds et alll2009l ). 

In this Letter we focus on SDSS J153636.22+044127.0 
(SDSS J153 6+0441 hereafter ) , a q uasar at a redshift 
z = 0.388 (|Boroson fc Lauerl I2009D . For the assumed 
flat cosm ology 3 , the inferred 0.1 p c separation subtends 
0.02 mas (Bo roson fc Lauerl l2009V This quasar has also 
been interpre t ed as a lone double-pe aked emitter (DPE; 
iGaskelll 120101: IChornock effall |2010h or as twin DPEs 



with a subparsec separation (jTang fc Grindlavl 2009). 
The emission lines from a DPE are thought to arise from 
rotational motion in a relativistic accretion disk. SDSS 
J1536 +0441 is radio quiet and our imaging of it at 8.5 
GHz ([Wrobel fc Laorl [20091 revealed two faint sources, 
SDSS J1536+0441A and J1536+0441B, separated by 
0.97" (5.1 kpc) with each source being unresolved with a 
diameter of less than 0.37" (1.9 kpc). It is now clear that 
SDSS J1536+0441A is hosted by the radio-quiet quasar 
(RQQ), while SDSS J1536+0441B i s hosted by a com- 
panion radio-loud ellip tical galaxy ([Decarli et al.l 120(391 : 
lLauer fc Borosonll2009D . 
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3 Ho = 71 km s — 1 Mpc - 1 and C m = 0.27, implying a luminosity 
distance of 2.1 Gpc, an angular size distance of 1.1 Gpc and a scale 
of 5.2 pc per mas. 



In the BBH scenarios for SDSS J1536 +0441A 
([Boroson fc Lauerll2009t pTang fc Grindlavil2009D . each of 
the two broad-line emission systems is itself a quasar. Re- 
cent optical spectroscopy finds that the broad-line emis- 
sion systems have relative positions that localize them to 
3 a bands of full wi dth 90 mas (470 pc) al ong a position 
angle (PA) of 48° ([Chornock et all 12010ft and 60 mas 
(310 pc) along PA 90° (T. Boroson, 2010, private com- 
munication). The resulting parallelogram localizes SDSS 
J1536+0441A to an area of 5400 mas 2 (0.15 kpc 2 ). Im- 
provements in the localization of the emission from SDSS 
J1536+0441A would further test its candidacy as a BBH 
system. In this Letter we use measurement techniques at 
radio frequencies to investigate the compactness of SDSS 
J1536+0441A, first by seeking evidence for a synchrotron 
self-absorbed spectrum and then by direct imaging with 
mas resolution to localize the emission on parsec scales. 
Our new imaging is presented in § [2] and its implications 
are explored in § [3] A summary appears in § 2] 

2. NEW IMAGING 

The CnB co n figura tion of the NRAO VLA 4 
([Thompson et al.l Il980f) was used under proposal 
code AL739 to observe SDSS J1536+0441A and 
J1536+0441B for 2 hours near transit on UT 2009 
May 27. We follo wed the strategies described by 
IWrobel fc Laorl ([20091 ) . except that for our new observa- 
tions the center frequency was 22.5 GHz, the switching 
time was 120 s, and the amplitude scale was set to an ac- 
curacy of about 5%. The elliptical Gaussian resolution, 
0.96" times 0.81" at PA -89°, was sufficient to obtain 
photometry for each source. We will report elsewhere on 
SDSS J1536+0441B. For SDSS J1536+0441A, the flux 
density was S = 1 .65 ±0.11 mJy and c omparison with 
the 8.5 GHz value ([Wrobel fc Laorll2009h implies a rising 

4 Operated by the National Radio Astronomy Observatory, 
which is a facility of the National Science Foundation, operated 
under cooperative agreement by Associated Universities, Inc. 
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spectrum with index a = 0.35 ± 0.08 (S oc v a ). The 
2009 December 31 release of the NRAO AIPS software 
was used for calibratio n and imagin g. 

The NRAO VLBA 4 (jNapier et al.lll994l) was used for 5 
hours under proposal code BL168 on 2009 October 14 UT 
to search for 8.4 GHz emission from SDSS J1536+0441A. 
Phase-referenced observations were made in the nodding 
style, using 32 MHz per circular polarization. A 120 s 
scan of SDSS J1536+0441A was preceded and followed 
by a 60 s scan of the reference source J1539+0430, fa- 
vorably located at a switching angle of 0.7° from SDSS 
J1536+0441A. The 2009 December 31 release of the 
NRAO AIPS software was used for calibration and imag- 
ing. We followed th e calibration strategies described by 
IWrobel fc Hoi (120061) and the la rge -area imaging strate- 
gies described bv lWrobel et alJ (|2005l ). 

We conducted a large-area VLBA searc h in a region set 
by th e 3 a VLA astrometric accuracy (|Wrobel fc Laorl 
120091) and of area 280,000 mas . Using natural weight- 
ing of the visibility data, this search achieved an ellip- 
tical Gaussian resolution of 2.57 mas times 1.04 mas 
at PA -5.2°, a geometric-mean resolution of 1.63 mas, 
and a root-mean-square sensitivity of 0.056 mJy beam -1 
in Stokes /. Given the large number of resolution el- 
ements sea rched, the con s ervat ive search strategies de- 
veloped by IWrobel et all (2005) were followed. Specif- 
ically, a 6 a detection threshold was adopted and only 
one source emerged above this threshold. This VLBA 
detection, shown in Figure 1, is unresolved with a di- 
ameter of less than 1.63 mas and has a flux density 
of S = 0.88 ± 0.12 mJy. The ratio of the flux den- 
sity o f the VLBA detection to that measured with the 
VLA (| Wrobel fc Laorl I2009D is 0.75±0.11. The position 
of the VLBA detection is a(J2000) = 15' l 36 m 36 i :2232 
and 5(J2000) = +04°41'27'/069, with a conservative 1 a 
error of 2 mas per coordinate estimated from analysis of 
the check source, J1544+0407, observed at a switching 
angle of 1.5° from the reference source. 

We also conducted a slightly deeper (5 a) small-area 
search for an additional point-like source near the VLBA 
detection and within the 5400 mas 2 localization region 
described in §[T]for the broad-line emission systems. Fig- 
ure 2 shows that if an additional source is present within 
that parallelogram-shaped region, it is fainter than the 
VLBA detection by a factor of 2.6 or more. This up- 
per limit was corrected for a 20% coherence loss, esti- 
mated from the ratio, about 1.2, of the VLBA detection's 
integrated-to-peak flux densities. This ratio is applica- 
ble to a switching angle of 0.7° and is also consistent 
with the higher ratio, about 1.6, measured for the check 
source which had a switching angle of 1.5°. 

3. IMPLICATIONS 

3.1. Size of the Radio-Sphere 

From the new VLA photometry, the RQQ SDSS 
J1536+0441A has a rising spectrum with index a = 
0.35 ± 0.08 at rest frequencies of tens of gigahertz. (This 
is consistent with our prior suggestion that the over- 
all spectrum of A and B was flat or rising between 1.4 
GHz [White et al. 1997] and 8.5 GHz [Wrobel & Laor 
2009].) SDSS J1536+0441A thus resembles the 45%-50% 
of RQQ that show flat or rising integrated spectra at 
similar frequencies (Barvain is et al.lll996l : lUlvestad et alJ 



120051) . This spectrum suggests that SDSS J1536+0441A 
is compact enough to be synchrotron self-a bsorbed, as ex- 
pecte d in the coronal framework for RQQ (|Laor fc Beharl 
l2008h . SDSS J1536+0441A has a bolometric luminos- 
ity o f L bol = 1.5 x 10 46 ergs s" 1 H oroson fc Lauerl 
2009) and a 22.5 GHz luminosity density of Lr — 
8.6 x 10 30 erg s s" 1 Hz" 1 . Applying equation (22) of 
lLaor fc Beharf (|2008t ). for a homo geneous synchrotron 
source with equipartition between magnetic and pho- 
ton energy densities, implies a radio-sphere of radius 
about 0.04 pc at a rest frequency of 31.2 GHz, and 
thus T b - 5 x 10 10 K. This is just below the iReadheadl 
(| 19941 ) limit of T b ~ 10 11 K, expected for equipartition 
between the electron energy density and the magnetic 
energy density within the radio-sphere. The 0.04 pc ra- 
dius corresponds to about 800 gravitational radii for a 
10 9 M BH, a mass thought to be applicable to SPSS 
J153 6+0441A (|Lauer fc Borosonll2009l ; iTang fc Grindlavl 
2009). The radius of the radio-sphere implies that the 
31.2 GHz emission arises just within the optical broad- 
line emission region (BLR) discussed in § 13.21 

The spectral index of this RQQ is clearly too shal- 
low to be a homogeneous source, and likely implies a su- 
perposition of emission from an inhomogeneous source, 
as commo nly adopted fo r flat spectrum radio-loud sys- 
tems (e.g.. lPhinnevlll985T ). Given the newly measured flat 
spectral i ndex for SPSS J1536+0441 A, its ratio of radio 
to X-ray ([Arzoumananian et al.|[2 009) luminosities drops 
from 5.9 x 10~ 5 (|Wrobel fc Laorl I2009D to 1.2 x 10~ 5 , 
putting it close to the 10 -5 average ratio characteriz- 
ing lower-luminosity active galactic nuclei, a ratio ex- 
pected in t he coronal framewor k for the radio emission 
from RQQ (|Laor fc Behanl2008l) 

Some RQQs are time variable ([Barvainis et al.f2005[ ) so 
our inference of a rising spectrum for SPSS J1536+0441A 
is weakened by using non simultaneous photometry. 
However, for the source to have an optically thin spec- 
tral index steeper than a = —0.5, the flux density at 
either 8.5 or 22.5 GHz would need to vary by a factor 
greater than two between the 100 days separating the 
measurements. Causality arguments would then imply a 
size upper limit of 100/(1+0.388) ~ 70 light days, or 0.06 
pc. In this case, both the steep spectrum and the high 
variability brightness temperature, T b > 10 10 K, would 
exclude a thermal free-free origin for the radio emission. 

From the new VLBA imaging (Figs. 1 and 2), a sin- 
gle source was detected at 8.4 GHz within the 470 pc 
by 310 pc localization region for the broad-line emis- 
sion (Chornock et al. 2010; T. Boroson, 2010, private 
communication). The VLBA detection has a geometric- 
mean diameter of less than 8.5 pc and a rest-frame 
brightness temperature, modified for an elliptical Gaus- 
sian, of T fc > 1.2 x 10 7 K. The ratio of the VLA and 
VLBA flux densities near 8 GHz is broadly consistent 
with unity. This suggests that the VLBA recovers all 
of the VLA signal but time variability remains a con- 
cern. For now, we tentatively assign the VLA-derived 
spectral index, a — 0.35 ± 0.08, to the VLBA detec- 
tion. Then the isotropic power at a rest frequency of 
8.4 GHz is P v = 2.9 x 10 23 W Hz" 1 . This VLBA 
detection has a power and T b limit at the low end of 
the val ues reported for othe r RQQs detected with the 
VLBA (iBlundell et all fl99& IBlundell fc Beaslevl Il99l 
lUlvestad et alJ 120051 ). The crude estimate made above 



Parsec-Scale Localization of SDSS J1536+0441A 



3 



for a synchrotron-self-absorbed size is also consistent 
with the VLBA detection. 

The implications of the above findings are examined 
below, first wit hin the context of BBH scenarios for SPSS 
J153 6+0441A (jBoroson fc Lauerll2009l : iTang fc Grindlavl 
I2009T) and then within the context of it being a lone PPE 
(|Gaskelll[2010t iChornock et al.l[2010h . 

3.2. Binary Black Hole Scenarios 

Our VLBA findings are consistent with the projected 
separation of the two quasars being less than 8.5 pc. The 
VLBA localization area for SDSS J1536+0441A is 82 pc 2 , 
improving over the emission-line localization area by a 
factor of about 1800. The VLBA detection appears to 
have a rising, synchrotron self-absorbed spectrum, which 
bodes well for imaging it at higher resolutions to improve 
the localization further. Optical spectroscopic moni- 
toring of SDSS J1536+0441A implies an orbital perio d 
longer than about 200 years (jLauer fc Borosorj 12009). 
Unfortunately, such a long period means that VLBA 
monitoring could not usefully constrain the astrometric 
wobble of SDSS J1536+0441A. 

The twin DPE model of (|Tang k Grindlavl I2009D re- 
produces the observed H/3 line profile by emission from 
a disk extending from 7000 down to 800 gravitational 
radii. Thus, the size of the radio-sphere at 31.2 GHz is 
just below the inner boundary of the BLR. Since the size 
of an optically thick radio-sphere scales as f„ A /v (e.g. 
equation 22 of Laor & Behar 2008), observations at mm 
wavelengths will allow us to probe the radio-sphere on 
smaller scales, potentially down to the optically emitting 
region at a few 10s of gravitational radii. If the source 
remains optically thick, then sub-mm observations can 
probe down to the X-ray emitting region at a few gravi- 
tational radii. 

The VLBA detection of the RQQ SDSS J1536+0441A 
represent the first parsec-scale localization of a can- 
didate BBH system identified through its broad-line 
properties. This VLBA detection also demonstrates 
that parsec-scale localizations of candidate BBH sys- 
tems need not be restricted to radio loud objects like 
the radio galaxy 040 2+379 with its 7-pc separation 
(jRodriguez et al.ll2006t ). 

Our VLBA findings cannot exclude additional sources 
with Tb < 4.5 x 10 6 K in the same field of view. 
Such a value is atypically low compared to other 
RQQs detected with the VLBA (IBlundell et all Il996t 
IBlundell k Beaslevl 119981: lUlvestad et al.l I2005D . But 
those studies targeted RQQ stronger than several mil- 
lijansky and were thus biased toward detecting higher 
brightness temp eratures. Several RQQ wer e not detected 
in the survey of IBlundell fc Beaslevl (fl998l ). with bright- 
ness temperature limits similar to the present study. This 
suggests that the effects of source resolution could also 
contribute to non detections of RQQs. 

3.3. Lone Double-Peaked Emitter Scenario 

The lone DPE scenari o for the RQQ S PSS 
J1536+0441A (|Gaskelll 120101: IChornock et~aHl201(l re- 
quires the presence of only one quasar. The line pro- 
files for SDSS J1536+0441A do make it an unusual 
DPE however (|Chornock et al.l 120101 : lLauer fc Borosonl 
2009). As a class, DPEs are rare, constituting only 4% 
of the spectroscopically-selected sample at z < 0.332 of 



IStrateva et all (|2003ft . Yet 76% of those DPEs are ra- 
dio quiet like SDSS J1536+044 1 A. Using 1.4 GHz de - 
tections from lWhite et al.l (jl997t ). IStrateva et ail (|2003f ) 
tabulate 1.4 GHz luminosities for 18 RQ DPEs and 
report typical values of several times 10 39 ergs s" 1 . 
SPSS J1536+0441AB is not detected by iWhite et al.l 
(1997) and the 1.4 GHz luminosity is less than 7.2 x 
10 39 ergs s -1 , a limit c onsiste nt with detected RQ 
PPEs in IStrateva efaLl (|2003t ). Moreover, except 
for SPS S J1536+0441A, no RQQs detected with the 
VLBA (IBlundell et all fl996l IBlundell fc Beaslevl H998I: 
lUlvestad et al.l l2005h are known to exhibit PPEs, so 
it is impossible to say whether or not the VLBA de- 
tection of SPSS J1536+0441A is in any way unusual. 
In this regard , VLB A imaging of the PPE sample of 
IStrateva et al.l (|2003| ) would be a useful undertaking. 
As the VLBA detection of SPSS J1536+0441A demon- 
strates, such imaging is feasible for both radio-quiet and 
radio-loud PPEs. 

4. SUMMARY 

Our VLBA search for 8.4 GHz emission from the RQQ 
SPSS J1536+0441A found only one source within the 
localization region, of area 0.15 kpc 2 , for the broad- 
line emission. The VLBA detection has a diameter of 
less than 8.5 pc and a T b > 1.2 x 10 7 K. This detec- 
tion of SPSS J1536+0441A represents the first parsec- 
scale localization of a candidate BBH system identified 
through its broad-line properties. The VLA photome- 
try at a rest frequency of 31.2 GHz yields a rising spec- 
trum, consistent with synchrotron self-absorption, which 
implies a radius of about 0.04 pc for the radio-sphere, 
and Xfc ~ 5 x 10 10 K. The observed compact flat spec- 
trum radio sphere is consistent with the trait predicted 
in the coronal framework for RQQs. The radio-sphere at 
31.2 GHz happens to be just inside the estimated inner 
boundary for the BLR, of 800 gravitational radii, in this 
object. 

It would be useful to investigate the spectrum at higher 
frequencies, corresponding to the mm and sub-mm range, 
to determine where the spectral slope steepens as the 
source becomes optically thin. This will establish the size 
of the most compact synchrotron-emitting region, and 
potentially allow a direct exploration of rclativistic elec- 
trons in the accretion disk corona. We plan such inves- 
tigations using the Expanded VLA ()Perlev et al.l l2009h 
and the Atacama Large Milli meter / submillimeter Array 
(|Wootten fc Thompsonll2009l ). 

Concerning the BBH scenarios for SPSS 
J1536+0441A, the VLBA detection is consistent 
with a quasar separation of less than 8.5 pc. No 
additional sources with Tj = 4.5 x 10 6 K or more are 
found in within the localization region for the broad-line 
emission. 

Concerning the lone PPE scenario for SPSS 
J1536+0441A, its emission line profiles do make it an un- 
usual PPE. But as no other radio-quiet PPEs have been 
imaged on parsec scales, it is impossible to say whether 
or not the VLBA detection of SPSS J1536+0441A is in 
any way unusual. 
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25 20 15 10 

MilliARC SEC 

Fig. 1. — VLBA image of Stokes I emission from SDSS J1536+0441A at a frequency of 8.4 GHz and spanning 20 mas (104 pc). The rms 
noise is 0.056 mjy beam -1 (1 <r) and the hatched ellipse shows the Gaussian beam dimensions at FWHM. Geometric-mean beam width is 
1.63 mas (8.5 pc) at FWHM. Contours are at —6, —4, —2, 2, 4, 6, 8, 10, 12, ... 14 times 1 a. Negative contours are dashed and positive 
ones are solid. Image peak is 0.79 mjy beam -1 . Linear gray scale spans —0.24 mjy beam -1 to 0.79 mjy beam -1 . 
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Fig. 2.— Stokes / emission at 8.4 GHz centered on the VLBA detection of SDSS J1536+0441A and spanning 200 mas (1.04 kpc). The 1 
a rms noise, beam dimensions, image peak and gray scale are as in Fig. 1. Contours arc at —5 and 5 times 1 a. The square shows the field 
of view for Fig. 1. The parallelogram encloses the localization region, of area 5400 mas 2 (0.15 kpc 2 ), for the broad-line emission systems. 
Any additional point-like source within the parallelogram has an observed peak below 5 a. 



